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1 Introduction 

Fractal percolation has been introduced by Mandelbrot in 1974 as a model 
for turbulence and is discussed in his book The Fractal Geometry of Nature 
f7\. Several equivalent formal definitions of this process can be found in 
the literature (see e.g. J2H2HH1). Here we only give an informal definition 
of the two-dimensional case. Let Kq be the unit square and choose an in- 
teger M > 2 and a parameter p G [0, 1]. To obtain K\, divide Kq into M 2 
equal subsquares, each of which survives with probability p and is dis- 
carded with probability 1 — p, independently of all other subsquares. Now 
do the same procedure in all surviving squares, in order to obtain A'2. Iter- 
ating this process gives a decreasing sequence of sets (K n ) n ^, see Figure 
[lj Let K = H neN K n be the limit set. 




Figure 1: Realizations of K n ,n = 1, 2, 3, 5, 7 for M = 3 and p = 0.85. 

It was shown in 1988 by Chayes, Chayes and Durrett [2] that there exists 
a non-trivial critical value p c (M) such that a.s. the largest connected com- 
ponent in K is a point for p < p c (M) and with positive probability there 
is a connected component intersecting opposite sides of the unit square for 

P > Pc(M). 

For all M > 2, the value of p c (M) is unknown. Several attempts have been 
made to find bounds for p c (M). It is easy to see that K is empty a.s. if 
P < 1/M 2 , which implies p c (M) > 1/M 2 . The argument in [2J is already a 
bit smarter: any left-right crossing has to cross the line {1/M} x [0, 1] some- 
where. A crossing of this line in K n means that there is a pair of adjacent 



squares on opposite sides of this line. Such pairs form a branching process 
with mean offspring p 2 M and consequently p c (M) > \j\f~M. For the case 
M = 2 this was sharpened by White in 2001 to p c (2) > 0.810, who used a 
set that dominates K and has a simpler structure to study. 
Sharp upper bounds are harder to obtain. The first idea to get rigorous up- 
per bounds for M > 2 was given by Chayes, Chayes and Durrett |2|, but (in 
their own words) these bounds are ridiculously close to 1. For M = 3, they 
show that p c (3) < 0.9999 (although in fact one can prove that p c (3) < 0.993 
with their method), which was improved by Dekking and Meester [4] to 
p c (3) < 0.991. Chayes et al only treat M = 3, but they point out that the 
same idea works for any M > 3. The case M = 2 can be treated by compar- 
ing with M = 4. As is noted by van der Wal a coupling argument gives 
Pc(2) < 1 — (1 — -y/p c (4) ) 4 . Following this approach gives p c (4) < 0.998 and 
p c (2) < 1 - 10- 12 . 

In this paper we present ideas to find significantly sharper lower and up- 
per bounds. To find lower bounds, we compare fractal percolation with site 
percolation. A fundamentally new result is that for all M we construct an 
increasing sequence that converges to p c {M). The terms in the sequence 
can in principle be calculated algorithmically. To actually obtain numerical 
results we bound them from below. These ideas lead to (computer aided) 
proofs thatp c (2) > 0.881 and p c (3) > 0.784. 

For the upper bounds, we map all possible realizations of K n to a finite al- 
phabet The choice of the alphabet is inspired by van der Wal's work 
19 j. The fractal percolation iteration process now induces an iterative ran- 
dom process on £/, which is easier to analyze than the original process. 
This theoretical framework is the basis of computer aided proofs for the 
following upper bounds: p c {2) < 0.993, p c (3) < 0.940 and p c (4) < 0.972. 

2 Lower bounds for p c {M) 

In this section we develop methods to calculate lower bounds for the criti- 
cal value of two-dimensional fractal percolation. First we briefly introduce 
site percolation and then we prove a coupling with fractal percolation that 
allows us to find lower bounds for p c {M). In particular, we construct an 
increasing sequence of lower bounds and we prove that this sequence con- 
verges to p c (M). At the end of this section we show how to use these in- 
sights to obtain numerical results. 

2.1 Site percolation 

Consider the infinite two-dimensional square lattice in which each vertex is 
open with probability p and closed otherwise. In this model the percolation 
probability ((p) is defined as the probability that the origin belongs to an 
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infinite open cluster. The critical probability is given by 



pf te :=inf{p:C(p)>0}. 

It has been shown by van den Berg and Ermakov (TJ that pf te > 0.556. 
Now consider the probability that a box of n x n vertices is crossed by an 
open cluster. This probability is monotone in n, so it converges to a limit. 
It is a classical result that a positive percolation probability is equivalent 
to this limit being positive (see e.g. [6| and the references therein). As a 
consequence we have the following property that will be used to couple 
site percolation to fractal percolation. 

Property 1 Take a box ofnxn vertices. Suppose p < pf te . Then the probability 
that there is an open cluster intersecting opposite sides of the box converges to as 
n — > oo. 

2.2 Coupling site percolation and fractal percolation 

In fractal percolation, a set in the unit square is said to percolate if it con- 
tains a connected component intersecting opposite sides of the square. Let 

6 n (p, M) = F(K n (p, M) percolates), 9(p, M) = F(K(p, M) percolates). 

Then p c (M) := inf {p : 8(p, M) > 0}. We will often suppress some of the 
dependence on M and p. It is well known that Hindoo 9 n (p) = 6(p) = 
P(lT=o i K n{p) percolates}), see 0. 

To obtain a proper coupling, we have to deal with diagonal connections. 
For example, the set [0, 1/2] 2 U [1/2, l] 2 percolates. We would like to ignore 
such diagonal connections, since in site percolation diagonal connections 
do not exist. Therefore we redefine percolation as follows. We say K n per- 
colates if there is a left-right crossing of the square that does not use di- 
agonal connections. Since diagonal connections break down almost surely, 
this does not change the limiting percolation probability. Also from now 
on connections, crossings and connected components in K n are similarly 
redefined. 

Theorem 1 Let ir n (p, M) = F(at least two sides are connected in K n (p, M)). If 
7r n (p) < pf te for some n, then p < p c (M). 

Proof First define delayed fractal percolation: F m>n is constructed in the same 
way as K m+n , the only difference being that we do not discard any squares 
in the first m construction steps. So we first divide the unit square into 
M m x M m subsquares and only then we start the fractal percolation process 
in each of these squares. Delayed fractal percolation stochastically domi- 
nates fractal percolation. 
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Figure 2: Three coupled realizations to illustrate ([TJ. Left: K$ for M = 2 and 
p = 2/3. Middle: F 2>3 for M = 2 and p = 2/3. Right: #i for M = 4 and 

p = vr 3 (2/3,2). 

Then we have the following inequalities (illustrated for M = m = 2 and 
n = 3 in Figure [2j: 

m+n (p, M) < P(F m ,„(p, M) percolates) < 0i(7r„(p), M m ). (1) 

The first inequality follows from the fact that i^ m , n is probabilistically larger 
than K m+n . The second inequality can be explained as follows. Suppose we 
have two types of squares: realizations of K n (p, M) in which at least two 
sides are connected (type 1) and realizations in which no side is connected 
to any other side (type 2). Suppose we tile a larger square with M 2m in- 
dependent realizations of K n (p, M). So the probability on a type 1 square 
is 7r n (p). This larger square is a (scaled) realization of -F m ,n(p, M). Now re- 
place all type 1 squares by a full square and discard all type 2 squares. This 
gives a realization of K\{tt(jp), M m ). Moreover, this replacement procedure 
can not destroy percolation. To see this, suppose we have a left-right cross- 
ing in F m>n (p, M). All parts of this crossing that are in type 1 squares are 
preserved. If the crossing uses a type 2 square S, then it enters and leaves S 
at the same side of S. At this side the neighboring square T must be of type 
1. So if we remove 5 and replace T by a full square, there still is a left-right 
crossing. 

A first level fractal percolation set can be seen as site percolation in a finite 
box. Suppose 7r n (p) < pf te for some n and let the box size M m tend to oo. 
By Property [T] we arrive at 

lim e 1 (7r n (p),M m ) = 0. (2) 

m— >oo 

Therefore, if 7r n (p) < pf te , by |lj we find that 

6(p) = lim 9 m+n (p) = 0, 

which in turn implies p < p c (M). □ 



4 



2.3 A convergent sequence of lower bounds for p c (M) 

In this section we define a sequence of lower bounds for p c (M). We prove 
that this sequence converges to the p c (M). Let 

p n c (M) = sup {p : 7T n (p, M) < pf e } . (3) 

Note that 7r n (p") < pf te for all n, since ir n (p) is continuous in p. Theo- 
rem [T] requires a strict inequality, but this does not give any problem since 
% n (p) is strictly decreasing in n, so 7r„ + i(p™) < p s c lte for all n and hence 
indeed p"(M) < p c (M) for all n. The strict monotonicity in n also im- 
plies that (p£(M))5£L is increasing. The obvious question now is whether 
(p"(M))^ =0 converges to p c (M). We will show that this is indeed the case. 
First we need that % n (p) goes to zero if the fractal percolation is subcritical: 

Lemma 1 Ifp < p c (M), then lim^oo 7r n (p, M) = 0. 

Proof Suppose p < p c (M), so 6{p) = 0. Note that a.s. there is an n such that 
in K n the two squares in the top left and bottom left corner are discarded 
already. Conditioned on this event, a connection in the limiting set from 
the left side to any other side can only occur if it horizontally crosses the 
vertical strip S consisting of the squares [0, M _n ] x [jM~ n , (j + 1)M _1 ] for 
j = 1, . . . , M n — 2. From self similarity and subcriticality it follows that in 
the limit each of these squares has zero probability to contain a component 
connecting opposite sides (horizontally and vertically). So in A" a horizon- 
tal crossing of S can only occur if it crosses a block of two vertically adjacent 
squares. But as Dekking and Meester showed (Lemma 5.1 in [4]), 9(p) = 
implies that such a block crossing has zero probability as well. It follows 
that connections from the left side to any other side have zero probability 
to occur. If this holds for the left side, then it holds for all sides, and so 
lim n ^oo7r n (p) =0. □ 

The previous lemma makes it easy to prove convergence: 

Theorem 2 The sequence (p™(M))^L converges top c {M) if ' n — > oo. 

Proof Let e > and suppose p = p c (M) — e. Then by Lemma [l] there exists 
an N such that 7r n (p) < pf te for all n > N. Therefore, (M) > p for all 
n > N. □ 

The theory developed so far gives us an algorithmic tool to calculate an 
increasing and converging sequence of lower bounds for p c (M). Still, it is 
not so easy to actually obtain sharp bounds. For example, a little thought 
gives 

7Tl(p,2) = l-(l-p) 4 , 

vri(p,3) = l + (l-p) 4 (p 5 + 4p 4 (l-p) + 6p 3 (l-p) 2 -l), 
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so using the bound of van den Berg and Ermakov for pf te , we find p\ (2) > 
0.183 and p\ (3) > 0.178. Of course, we want to find a bit sharper bounds, so 
we should take larger values of n. However, for large n, the functions TT n (p) 
are very complicated polynomials, and it is not clear how to find them in 
reasonable time. In the next section we will discuss a way to avoid this 
problem. 

3 Classifying realizations 

In this section we introduce the idea to map realizations of K n to a finite 
alphabet si that does not depend on n. Such a map will be called a classifi- 
cation. In this way, we can simplify the substitution process without losing 
too much essential information on the connectivity structure in K n . The 
construction of the random sets (K n )^L now induces a sequence of prob- 
ability measures on si '. These ideas can be used to obtain both lower and 
upper bounds for p c (M). For the lower bounds we choose the classifica- 
tion in such a way that we can calculate upper estimates for n n (p) and then 
we can approximate terms of (p™)™ =0 using the van den Berg and Ermakov 
bound for p^. lte . For the upper bounds we will construct a classification that 
permits us to bound 9{p) away from for p sufficiently large. 
Let J(f n be the set of all possible realizations of K n . We will define a par- 
tially ordered alphabet si and the classification will actually be a sequence 
of maps c € = ( < ^„)^ =0 , where ^„ is a map from J^ n to si. The alphabet will 
have a unique minimum and maximum, with the property that 

^f„(0) = min(^) and tf n ([0, l] 2 ) = max(j^). (4) 

Realizations of K n can be mapped in a natural way to words in si . 
A realization of K n is obtained by first generating a realization of K\ and 
then replacing all squares that survived by independent scaled realizations 
of K n -\. Given ^ n -i, this induces a map W n from Jif n to si MxM . If ^„(K n ) 
only depends on the word W n (K n ) and not on n, we say the classification 
is regular. 

Definition 1 Let <pbea map from si MxM to si. Define = ( < ^ 7 n )^ =0 by letting 

tf n = ( f ) oW n , n>l. (5) 

Then we say is a regular classification and <f> is called the word code ofr^. 

Example 1 Let M = 2 and take the alphabet si = {□, ffl}, where min s/ = □ 
and max s/ = EEL For a 2 x 2 word w over si let 4>(w) = EB if and only if at 
least one of the letters in w equals EEL Let ^ = ( < ^ ri )^ be regular with word 
code (p. Then % is determined by Q: 

%(K o (p)) = %{[0, l] 2 ) =EB. 
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For n > 1, = <j) o W n . For instance 

«b(0) %([o,i] 2 ) i _ , r □ a 1 

«b(0) «b(0) j 1 □ □ / " 



We now want to analyze the probabilities P^ (-fT n (p)) = a), where a G jz/. 
Suppose ^ is a regular classification with word code 4>. Let x G [0, l]'^' be 
a probability vector on srf, and suppose we construct an M x M word u; in 
which all letters are chosen independently according to x. Define F^(x) G 
[0,l]Klby 

(i f V(a;))o = Px(0(w)=a), se/ 

The function F v : [0, l]^ -)• [0, l]^ will be the key to calculate the proba- 
bilities P(^ \K n {p)) = a) in an iterative way as is shown in the next lemma. 
Let r n (p) be a vector indexed by the elements of si such that 

r:(p)=F(^(K n (p)) = a). 

Let r D and r ffl be the vectors that assign full probability to min(^) and 
max(j/) respectively. 

Lemma 2 iff foe classification is regular, then 

r n+1 (p) = Ftf(pT n (p) + (1 — p)r a ) with initial condition r°(p) = r ffl . 

Proof The letters in W n+ \{K n+ i(p)) are independent. With probability p a 
letter corresponds to a scaled realization of K n (p), with probability 1 — p 
it corresponds to an empty square. So each letter occurs according to the 
probability vector pr n {p) + (1 — p)r a . This gives the recursion. The initial 
condition follows from Q. □ 

This recursion formula is essentially a generalization of the recursion given 
inH. 

3.1 Strategy for lower bounds 

The strategy to find lower bounds for p c now is as follows. Define an al- 
phabet ff/ with subset &Z W and define a classification ^ (by choosing 0) in 
such a way that F(^f(K n (p)) G jz4) > ^n(p) f° r all n - Now search for the 
largest p for which this probability is smaller than 0.556. We will give an 
example to illustrate this procedure. In this example (which only gives a 
very moderate bound) only has two elements, and therefore exactly the 
same recursion as in [4j pops up. 
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Example 2 (A lower bound for M = 2) Let ^ be the classification of Ex- 
ample [I] By induction it follows that if two sides are connected in K n , then 
^ n (K n ) = EEL Consequently r^(p) > 7r n (p). From the definition of it fol- 
lows that in this case 

r a " +1 (p) = 1 - (1 -pr«(p)) 4 = ^(pr^p)). 

The function G p (x) : = (pa;) is increasing and (p) decreases to the largest 
fixed point of G p (-). Choosing p = 0.33, we find tt^,q{p) < t|°(p) m 0.554 < 
0.556 < p s c lte , and consequently pf(2) > 0.33 and therefore also p c (2) > 
0.33. ■ 



3.2 Strategy for upper bounds 

Our recipe to find upper bounds for p c {M) is a bit more involved. We al- 
ready have a partial ordering on srf , but to find upper bounds it will also 
be convenient to partially order the probability vectors on srf . A set S C $rf 
will be called increasing if a G S implies b G S for all b y a. For proba- 
bility vectors x and y, we now write x >z y if x assigns larger probabili- 
ties to all increasing subsets of s& '. We say the function F<# is increasing if 
F<g{x) y Fc g {y) for x y y. 

Lemma 3 Let ^ be a regular classification for which Fcg is increasing. Then 
(r n (p))^ =0 is decreasing and r°°(p) := Hindoo r n (p) exists. If F^(px + (1 — 
p)t d ) y xfor some probability vector x G [0, l]^ and p G (0, 1], then r°°(p) y 
x. 

Proof Since r°(p) = r ffl , we have r°(p) ^ r 1 (p). Suppose r n (p) ^ r n+1 (p). 
Then r n+1 (p) = Fg (pr n (p) + ( 1 —p) t d ) ^ F^(pr n+1 (p) + (l-p)r D ) = r n+2 (p) 
since is increasing. So (r n (p))^L is decreasing. If S C ^ is an increasing 
set, then 

r£(p) := P(^(A' n ) G 5) 

decreases in n and is bounded from below by 0, so lim^oo r^(p) exists. 
There are only finitely many increasing subsets of srf and their limiting 
probabilities uniquely determine r°°(p). 

For the second statement, r°(p) = r ffl y x. Now suppose r n (p) y x for 
some n. Then 

r n+1 (p) = F v (p T n {p) + (1 -p)r D ) ^ F^(px + (1 -p)r D ) b z. 

and therefore r n (p) ^ x for all n G N. Hence T°°(p) ^ x. □ 

Corollary 1 Let "^bea regular classification for which Fcg is increasing. Suppose 
x G [0, l]!^ is a probability vector for which x mSbx{s/) > 0. 1/(9 n (p) > r^^O) 
/or flZ/ n, and F^(px + (1 — p)r D ) ^ x/or some p G (0, 1], ffaen 

0(p) = lim n (p) > lim r max( ^ } (p) > x max( ^) > 0, 
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and consequently p > p c . 

Before we give an example of the procedure to find upper bounds for p c (M), 
we will first construct suitable alphabets in the next section. 



4 Construction of the alphabet and word codes 



Our alphabets will be defined by means of non- 
crossing equivalence relations. Figure [3] displays 
all equivalence relations on an ordered set of 5 
elements, of which 42 are non-crossing. Let E = 
{ei, . . . , e„} be a collection of line segments that 
together form the boundary of the unit square, 
numbered clockwise, starting from (0, 0). We de- 
fine the alphabet as the set of all non-crossing 
equivalence relations on E. The size of s^e is given 
by the Catalan number ^ ( 2 ™) . 
On the alphabet s$e, we can naturally define a 
partial ordering. For a, b G s#e we say a >z b if 
all equivalences that hold in b also hold in a. One 
should think of the equivalences (denoted by ~) 
as connections between the elements. The letters 
can be represented as pictures in a square. For ex- 
ample, if E is the set containing the four sides of 
the square (so n = 4), then s^e can be represented 
as 

r . Figure 3: Equivalence 

b i'''''''''''''/ relations on 5 elements 

The alphabet has a unique maximum and minimum, max(^4;) = EH and 
m.m.(s$E) = From now on we will denote the alphabet that corresponds 
to dividing each side in M n elements of equal size by 
Now that we have letters, we can start making words. We only allow square 
words, obtained by tiling [0, M] 2 by M 2 letters. An example of a 3 x 3 word 
over £^3 o is given in the right panel of Figure [4] A word w also has a bound- 
ary element set E w . We enumerate clockwise, starting from (0, 0), as fol- 
lows: E w = {ef,...,e^ lE] y 
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4.1 Weak and strong connectivity 

Suppose we partition E w into subsets E^, . . . , E* for some integer k > 2. 
If e ~ / for some e G E l w and / G we say E l w and are connected. If 
£^ , . . . , E%p form a chain of pairwise connected sets, we say E 1 ^ and E l ™ 
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Figure 4: Left: A realization of K2 for M = 3. Right: A 3 x 3 word w over 
s^zfl- If we use the word code ^3,0, then this word corresponds to the real- 
ization on the left. 

are weakly connected. If strictly more than half of the elements in E l w and 
strictly more than half of E 3 W are in the same equivalence class, we say E l w 
and E ] w are strongly connected. Observe that weak and strong connectivity 
are equivalence relations on {E^, . . . , E^\. 

Choose the alphabet s^M,m then \E\ = 4M n , so words of size M x M have 
4M n+1 boundary elements. Partition E w in the subsets 

K = {eZ (i -i )+j :l<j<M], l<i<AM n . 

This gives a tool to define word codes: a word defines an equivalence re- 
lation on { E^ , E^ Mn } that obviously can be mapped to an equivalence 
relation on E = {e±, . . . , e^u™}, which is just a letter in s^e- In this way 
define word codes 

$M, n :^M,n M ^^M,n and * M) n : ^M,n M ~> ^M,n 

based on weak and strong connectivity respectively. The corresponding 
classifications are increasing. 

Example 3 Let w be the word in Figure |4j This word has 12 boundary el- 
ements that will be grouped into four partition sets. For instance, con- 
tains the three elements at the right side. Then E^, E^ and E^ are all pair- 
wise connected and E^ is connected to E^. Consequently E l w is weakly 
connected to E 3 W for all i and j. Therefore 3>3,o(w) = EE 
The only partition sets that are strongly connected are E^ and E^. There- 
fore ^>3fl(w) = □. ■ 
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The idea behind the definitions of <3? M,n and ^M,n is that they guarantee the 
following key properties: 

Property 2 Using the alphabet s^M,k and word code <&M,k> the following implica- 
tion holds: ifei, ej G E are connected in K n , then e% ~ ej in &(K n ). 

Property 3 Using the alphabet s$M,k an d word code ^m,u> the reversed implica- 
tion holds: ifei ~ &j in ^(Kn), then ej and ej are connected in K n . 

These properties are not hard to prove. Therefore we do not give a fully de- 
tailed proof. For Property [2] a straightforward inductive argument suffices. 
For Property [3j one should note that we actually have a stronger impli- 
cation, namely that and ej are connected in a special way that we will 
explain using the realization in Figure [4] In this realization the left side 
and the top side are connected. At both these sides there are two (strictly 
more than M/2) first level squares each of which contains two second level 
squares that survived. These second level squares are all in the same con- 
nected component. This property is easily generalized to higher level re- 
alizations and other values of M. If two neighbouring squares both have 
such connection from their joint side to an other side, then these two con- 
nections are in the same connected component. 

The following lemma shows that these alphabets and word codes are suit- 
able for our purposes, see the discussion in Section [3] 

Lemma 4 Using the word codes <&u,k an d ^M,k leads to two inequalities: 

1. Take the alphabet s#M,k an d define a classification *rf by the word code $ M,k- 
Let =«4 ^ s^M.h be the set of all letters in which at least two of the sides are 
connected. Then r"(p) := P(<;f (-fT n (p)) G ^4) > n n (p) for all n. 

2. Take the alphabet s^M,k and define a classification ^ by the word code ^M,k- 
Then r maxK)(p) ^ d n(p)for all n. 

Proof These statements follow from Property [2] and [^respectively □ 

Now we are ready to give an example illustrating how to find upper bounds. 
We will keep the example as simple as possible, so that it can be checked 
by hand. Therefore our alphabet will contain only two letters and we will 
use a word code that is even simpler than ^3,0- Nevertheless, it leads to a 
bound that already improves upon the best bound known so far. 

Example 4 (An upper bound for M = 3) Let srf = {max(^^), mir^^o)} = 
{□, ffl}. A 3 x 3 word w has twelve boundary elements that we partition into 
four sets as before. Define a classification ^ by letting 4>(w) = ffl if and only 
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if all sets E^, . . . , E^ are strongly connected. This classification is increas- 
ing and satisfies Property |3l so t^w^ (p) < On (p) for all n. The recursion 
of Lemma [^reduces to (p) = r^(prg(p)), where 

7ffl(p) = P 9 + 9p 8 (l - p) + 20p 7 (l - p) 2 . 

This function is the first component of F^(pr m + (1 — p)r D ) = i*V((p, \ —p)). 
Now choose p = 0.984 and x = 0.9720. Then r^(pz) ^ 0.9721 > x. There- 
fore the probability vector (x, 1 — x) satisfies the requirement of Corollary 
[l] Hence p c (3) < 0.984. ■ 

4.2 Monotonicity and convergence 

So far we developed some tools to find bounds for p c (M). One would ex- 
pect that taking larger alphabets results in sharper bounds, since we can 
approximate the connectivity structure in K n more accurately. In this sec- 
tion we show that this is indeed the case and that the lower bounds even 
convergence to p c {M). Unfortunately, we do not know if the upper bounds 
also converge to p c (M). 

For the word code &M,k over s^M,k, define the corresponding classification 
and let r"(p) be defined as before. Then define a critical value as follows: 

p c {®M,k) :=sup{p:C(p)<pf e }. 

Let srf = s^M,k and define ^ by the word code ^M,k- Also here we define a 
critical value: 

Pc(^M,k) := mf {p : C^p) > 0} . 
Now we have the following proposition: 

Proposition 1 The sequence (p c ($M,k))kLo zs increasing and (p c (* 'M,k))kLo zs 
decreasing. Moreover, 

lim p c ($M,fc) =Pc(M). 

k— >oo 

Proof The alphabet k contains equivalence relations on an element set 
Eu,k- The element set EM,k+l of £^M,k+i can be obtained by dividing each 
element e € Em,h m t° M equal pieces. These pieces will be called children 
of e, and e will be called the parent of its children. Define classifications 
^ k and by the word codes §M,k an d $M,k+l- Denote the correspond- 
ing probability vectors by k r and fc+1 r. By induction it follows that if two 
elements are connected in < ^ 7fc+1 (iir n ), then their parents are connected in 
^ k (K n ). Consequently, fc r"(p) > k+1 T™(p) and so 

Pc (*(M,fc + l)) >p c ($(M,fc)). 
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Now define ^ and V k+1 by the word codes *M,fc and ^M,k+i- If f wo e l e ~ 
ments are connected in < tf k (K n ), then two of their children are connected in 
^ k+1 {K n ). Henceforth, fc r« ax ,» < k+1 r^(p), so 

Pc (*(M,k + l)) <p c {V{M,k)). 

Take the alphabet s#M,k an d define ^ k by the word code <1> A realization 
of K n consists of M n x M n squares, and letters in £^M,k have M k boundary 
elements at each side. This means that for n < k the classification describes 
the connectivity structure exactly: two elements in (K n ) are connected 
if and only if they are connected in K n . So ir n (p) = fc 7""(p) if n < k. This 
implies that for n < k we can rewrite Q: 

p n c (M) = sup {p: k r^p)<pf e } 

< sup [p : *r~(p) < pf e ) = p c ($ M ,k) < Pc(M), 

where we used that k r™(p) decreases in n by Lemma [3] Theorem [2] states 
that Pc(M) converges to p c (M), so we conclude that lim^oo Pc(&M,k) = 
Pc(M). ' □ 



5 Numerical results 

In this section we present our numerical results. Our principal goal was 
to compute bounds using the alphabets £^u,k which were constructed in 
the previous section. In these calculations we encountered the problem of 
accumulating rounding errors, disturbing the convergence. Therefore we 
had to normalize the probability vector in each step. We did the rounding 
and normalization in such a way that our conclusions are not violated. For 
example to show that T™(p) < pf te , we made sure that we computed an up- 
per estimate for r™(p). Our implementation in Matlab (everything available 
from the author on request) gives the following results: 

Proposition 2 Take the alphabet s$M,k an d define a classification *rf by the word 
code <3?A/,fc- Let n = 1000 and define r™(p) as before. Then 



- For M = 


2 and k = 0, w e have r"(0. 


.785) 


^ Pc 


- For M = 


2 and k = 1, we have r™(0. 


.859) 


^ ~-.site 
~ i J c 


- For M = 


3 and k = 0,we have r™(0. 


.715) 


^ Pc 


Corollary 2 p c 


: (2) > 0.859 and p c (3) > 


.715. 





13 



20 40 60 80 100 

Figure 5: Plot of r"(p) for p = 0.7 + 0.01/c where fc = 0, . . . , 9 as functions 
of n. Note that there is a kind of phase transition between p = 0.78 and 
p = 0.79. 

Proof This follows from Lemma H] and Theorem [T] □ 

Figure [5] illustrates for the case M = 2 and A; = how r" (p) behaves as 
a function of n for some values of p. The values of r" (p) were calculated by 
iterating the recursion of Lemma |2j 

For larger values of k the computations were too complicated to perform in 
a reasonable computation time. For example, the alphabet 3/2,2 already con- 
tains 35357670 letters. Nevertheless we will show that it is possible to im- 
prove the bounds of Corollary |2]by taking other alphabets or word codes. 
For M = 2, define the element set E by dividing the left and right side of 
the unit square into 4 equal pieces and the bottom and top side into 2 equal 
pieces. This leads to an alphabet s/e that is in some sense in between 3/2,1 
and 3/2,2- Analogous to our previous approach, we define a classification 
by choosing a word code based on weak connectivity. For this classification 
we find (0.876) < pf e , which implies p c (2) > 0.876. 
One can improve this even a bit more by taking the alphabet 3/2,2 and defin- 
ing a word code $2,2 that is a bit simpler than $2,2 as follows. If at least one 
of the letters in a 2 x 2 word w equals 111111(^2,2)/ then &2,2(w) = &2,2{w). 
Otherwise, define $2,2 (w) by first mapping each of the four letters to s/e 
and then mapping the new word to 3/2,2, in both steps using weak connec- 
tivity. This simplifies the required calculations a lot, and leads to p c (2) > 
0.881, in 200 iterations. 
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For M = 3 we improved the lower bound by using an alphabet that divides 
the left and right side of the unit square into 3 equal pieces. The bottom and 
top side still consist of one element. This leads to p c (3) > 0.784, using 100 
iterations. 

These calculations have been checked by Arthur Bik, a mathematics stu- 
dent at Delft University of Technology. He independently implemented the 
algorithms and reproduced all results, except the bound p c (2) > 0.881. This 
was due to the fact that his program was not fast enough to perform the 
calculations in a reasonable time. Concluding, the best lower bounds we 
found are 

Theorem 3 p c (2) > 0.881 and p c (3) > 0.784. 

Now let us turn to the upper bounds. Here we want to make use of Corol- 
lary [Tl We already know ''^/^(p) < n (p) for all n (Lemma |4l and that 
the classification is increasing. That means, we only have to find a proba- 
bility vector x G [0, l]^ for which x maK r^\ > and p G (0, 1] such that 
F^(px + (1 — p)r a ) >z x. The trick we use here is to take a value of p, and to 
search for the fixed point r°° by iterating the recursion of Lemma|2] Assume 
that the numerical results suggest that T^ ax ^ > 0. Letting x = r n (p) for n 
large, we have a combination of x and p for which F<#(px + (1 — p)r D ) « x. 
Taking p slightly larger gives a combination for which the desired inequal- 
ity Ftf(px + (l— p)r a ) y x holds. This strategy leads to the following results: 

Proposition 3 Take the alphabet s^m,u an d define a classification by the word 
code ^M,k- Let n = 1000. The conditions F^(px + (1 — p)r D ) y x and x max r^\ > 
hold ifx and p are chosen as follows: 

- For M = 3 and k = 0, choose p = 0.958 and x = r n (0.9579). 

- For M = 4 and k = 0, choose p = 0.972 and x = r"(0.9719). 

Corollary 3 p c (3) < 0.958 and p c (4) < 0.972. 

Proof This follows from Lemma |4]and Corollary [T] □ 

For M = 3, the result can be sharpened by using the same alphabet that was 
used to improve the lower bound. The classification is again defined by a 
word code based on strong connectivity. In that case the choice p = 0.940 
and x = r 1000 (0.9399) satisfies all conditions, so p c (3) < 0.940. 
The algorithm for M = 4 can be slightly adapted to find a bound for M = 2. 
Each realization of K n for M = 4 can be seen as a realization of Km for 
M = 2. So the only thing that changes is the iteration function F<g. Given 
a probability vector on the probability on each 4x4 word can be cal- 
culated. The word code is still the same function ^4,0- These ingredients 
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are sufficient to determine F<g. Performing the calculations we find that the 
conditions are satisfied for p = 0.993 and x = r 1000 (0.9929) / henceforth 
p c (2) < 0.993. 

Also for the upper bounds Arthur Bik checked our results. He indepen- 
dently reproduced our bounds, except for the bound p c (3) < 0.940 (for 
similar reasons as before). Summarizing, our best upper bounds are 

Theorem 4 p c (2) < 0.993, p c (3) < 0.940 and p c (4) < 0.972. 
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